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Abstract 

A  simulation  model  for  a  tubular  solid  oxide  fuel  cell  (SOFC)  was  developed  by  the  object-based  approach  to  calculate  the  current 
distribution,  gas  concentration  distribution,  and  temperature  distribution  at  the  steady  states  and  transient  operation  states.  The  transient 
electrical  and  temperature  response  to  a  load  change  was  simulated  for  the  both  cells  with  the  diameter  of  22  mm  (standard  cell)  and  the 
diameter  of  2.4  mm  (micro-tube  cell).  The  time  required  to  reach  the  new  steady  state  as  the  operating  voltage  was  changed  from  0.7  to  0.5  V 
for  the  micro-tube  cell  was  found  to  be  15  s,  which  is  much  shorter  than  that  of  the  standard  cell. 
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1.  Introduction 

Solid  oxide  fuel  cell  (SOFC)  has  many  advantages  due  to 
its  high  operation  temperature,  such  as  high  energy  conver¬ 
sion  efficiency,  flexibility  of  usable  fuel  type,  and  high 
temperature  exhaust  gas.  However,  the  rapid  start-up  of 
SOFC  is  difficult  due  to  its  high  operation  temperature, 
so  a  good  promise  of  the  SOFC  system  as  a  stationary  power 
supply  system  was  mainly  surveyed  [1,2].  Recently,  it  is 
proved  that  SOFC  with  micro-tube  cell  having  diameters  of 
2  mm  could  be  heated  up  to  operation  temperature  within  a 
few  seconds  without  any  crack  [3].  Further,  micro-tube  cells 
can  reduce  the  system  volume  and  therefore  increase  the 
output  power  density,  because  the  effective  electrode  area 
per  unit  volume  increases  by  decreasing  the  cell  diameter. 
The  higher  output  power  density,  in  conjunction  with  the 
high  tolerance  to  thermal  stress  makes  it  possible  to  apply 
SOFC  systems  as  the  power  supply  to  vehicles  [4,5], 

During  the  operation  of  an  SOFC  system,  the  suppression 
of  mechanical  stress  and  the  quick  system  response  to  a 
demand  change  are  important.  Mechanical  stress,  which 
leads  to  the  mechanical  breakdown  of  the  cell  components, 
can  be  classified  into  two  types,  i.e.  residual  thermal  stress 
through  fabrication  and  thermal  stress  caused  by  the  tern- 
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perature  distribution  along  the  cell.  In  order  to  design  the  cell 
and  stack,  which  can  prevent  the  mechanical  breakdown,  the 
analysis  of  the  temperature  distribution  is  first  necessary. 
Further,  to  apply  SOFC  systems  to  vehicles,  it  is  important  to 
investigate  the  transient  characteristics  of  the  SOFC. 

Although  such  analyses  by  experimental  approach  are 
necessary,  it  is  time  consuming  to  investigate  the  character¬ 
istics  for  all  the  cells  at  various  operating  conditions  only  by 
experimental  approach  to  optimize  the  system  design.  Com¬ 
putational  modeling  is  an  effective  tool  for  designing  the 
optimal  SOFC  system  and  for  predicting  their  characteristics 
at  the  steady  states  and  transient  operating  states.  The 
transient  electrical  response  of  a  planar  cell  [6]  and  a  tubular 
cell  [7,8]  has  been  studied  using  the  conventional  cell 
component  materials.  However,  we  are  not  aware  of  any 
previous  literature  reporting  the  transient  characteristics  of 
SOFC  for  different  cell  and  stack  configurations  and  cell 
component  materials  on  the  same  basis. 

While  there  is  a  strong  need  to  compare  SOFC  with 
different  cell  configuration  and  different  cell  components, 
an  analysis  and  designing  of  the  SOFC  system  requires 
interdisciplinary  collaborations.  An  object-based  modeling 
approach  can  work  effectively  for  the  interdisciplinary 
collaboration  [9].  In  an  object-based  model  framework,  each 
component  model  is  modularized  as  an  object  with  explicit 
interface.  By  defining  the  interfaces  explicitly,  each  modeler 
only  needs  to  assure  that  the  interface  of  his  or  her  own 
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Nomenclature 

Lists  of  symbols 

C  concentration  (mol  m-3) 

Cp  heat  capacity:  solid  (J  kg  1  K-1); 

gas  (J  mol-1  K-1) 

D  diffusion  coefficient  (m2  s-1) 

E  activation  energy  (J  mol-1) 

emf  electromotive  force  (V) 

/  view  factor  (-) 

F  Faraday  constant  (=96485  C  mol-1) 

h  heat  transfer  coefficient  (W  m-2  s-1) 

Hm(T)  enthalpy  of  gas  species  m  at  T  (J  mol-1) 
H®n(T)  heat  formation  of  gas  species  m  at  T 
(J  mol-1) 

i  current  density  (A  m-2) 

Jm  molar  flow  rate  of  constituent  m  (mol  s-1) 

k  thermal  conductivity  (W  m-1  K-1) 

M  molecular  mass  (kg  mol1) 

p  pressure  (Pa) 

p0  total  pressure  (Pa) 

q  heat  (W) 

r  mean  pore  radius  (m) 

R  gas  constant  (=8.314  J  mol-1  K-1) 

S  surface  area  (m2) 

t  time  (s) 

T  temperature  (K) 

V  volume  (m3) 

Greek  symbols 

dx  width  of  each  sliced  cell  (m) 

e  porosity  (-) 

y  proportionality  constant  for  exchange  current 

density  (A  m-2) 

;;  overpotential  (V) 

X  emissivity  (— ) 

p  density  (kg  m-3) 

cr  Stefan-Boltzmann  constant  (W  m-2  K-4) 

£  tortuosity  (-) 

,9  transfer  coefficient  (— ) 

Indices 

a  anode 

act  activation 

at  air  feed  tube 

c  cathode 

cone  concentration 

cond  conduction 

conv  convection 

e  electrolyte 

fuel  fuel 

g  gas  (fuel,  oxidant) 

gen  heat  generation  by  entropy  change 

H2  hydrogen 


H20  water  vapor 

ox  oxidant 

02  oxygen 

rad  radiation 

react  enthalpy  gain  due  to  mass  transfer 

s  solid  (cell,  air  feed  tube) 

s-g  interface  between  solid  and  gas 


model  matches  the  interface  of  the  framework  instead  of 
understanding  all  the  implementation  details  of  other  mod¬ 
els.  This  enhances  the  flexibility  of  model  extension  and 
makes  it  easy  to  replace  models,  which  compose  the  model 
framework,  with  other  models  such  as  more  precise  models 
or  models  that  is  based  on  the  different  assumptions.  As  a 
first  step  to  model  SOFC  system  in  an  object-based 
approach,  we  have  developed  a  simulation  model  that 
simulates  the  characteristics  of  a  single  tubular  SOFC.  By 
using  the  developed  model,  time  responses  of  cells  and 
transient  temperature  distribution  profiles  for  both  the  stan¬ 
dard  tubular  cell  and  micro-tube  cell  are  compared  in  this 
study. 

2.  Model  description 

2.1.  Cell  configuration 

Fig.  1  illustrates  the  configuration  of  the  cathode-sup- 
ported  type  tubular  cell.  Zr().84Yo.i602_s  (YSZ),  Ni-YSZ 
cermet,  La0.81Sr0.09MnO3_s  (LSM),  and  La0.9Sr0.iCrO3_8 
(LSCr)  was  used  for  electrolyte,  anode,  cathode,  and  inter¬ 
connect,  respectively.  Although  it  is  not  explicitly  shown  in 
Fig.  1,  Ni  felt  was  used  for  the  electrical  contact  between 
cells.  The  physical  properties  of  each  component  material 
and  gas  are  shown  in  Tables  1  and  2  [10-13].  We  assumed 
that  resistivites  follow  the  Arrhenius  equation.  Because  the 
molar  specific  heat  of  LSM  was  not  found,  we  set  the  value 
as  100  J  mol-1  K-1  in  the  calculation.  In  this  study,  we 
evaluated  characteristics  of  SOFCs  with  two  different 
dimensions.  One  is  the  cell  with  the  diameter  of  22  mm 


Fig.  1 .  Configuration  of  tubular  SOFC. 
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Table  1 

Physical  properties  of  solid 


Component 


Anode 


Electrolyte 


Cathode 


Interconnect 


Air  feed  tube 


Material 

Ni-YSZ 

YSZ 

LSM 

LSCr 

A1203 

Resistivity  at 

8.14  x  1CT4 

45 

0.014 

3 

- 

1073  K  (Q.  cm) 

Activation  energy 

11.6 

-76.2 

-5.0 

0.12 

for  resistivity  (kJ  mol-1) 
Density  (g  cm-  ) 

6.87 

5.90 

6.57 

3.95 

Molar  specific 

52 

73 

100 

- 

120 

heat  (J  mol-1  K-1) 

Thermal  conductivity 

4 

4 

4 

10.4 

(WnP1  K~‘) 

Table  2 

Physical  properties  of  gas 


Component 

Molar  specific  heat 
(J  moP1  K_1) 

Thermal  conductivity 
(W  rrP1  K_1) 

Fuel 

- 

0.46 

Oxidant 

- 

0.071 

h2 

31 

- 

h2o 

43 

- 

02 

35 

- 

N2 

33 

- 

(standard  cell)  [14],  and  the  other  is  the  cell  with  the 
diameter  of  2.4  mm  (micro-tube  cell)  [4],  The  dimensions 
of  each  cell  and  each  cell  component  are  listed  in  Table  3. 

2.2.  Oveniew  of  the  model  framework 

Numbers  of  phenomena,  such  as  reaction,  mass  transfer 
and  heat  transfer,  and  physical  properties  influence  the 
characteristics  of  SOFC  system.  We  make  many  simplifica¬ 
tions  when  we  construct  simulation  models,  because  it  is  too 
unrealistic  to  construct  a  model  that  takes  all  the  phenomena 
occurring  in  the  system  into  consideration.  In  addition,  some 
of  the  physical  property  data  are  sometimes  not  available. 
Therefore,  it  is  important  to  develop  a  model  framework  that 
enables  users  to  flexibly  change  assumptions,  physical 
properties,  and  cell  configuration.  As  a  first  step  to  develop 
a  model  framework  that  satisfies  these  requirements,  we 
have  developed  a  model,  which  simulates  current  and  tem- 


Table  3 

Cell  dimension  for  the  standard  cell  and  the  micro-tube  cell 


Standard  cell 

Micro-tube  cell 

Effective  cell  length 

150  cm 

25  cm 

Outer  diameter  of  the  cell 

22  mm 

2.4  mm 

Thickness  of  cathode 

2  mm 

130  pm 

Thickness  of  anode 

100  |rm 

50  pm 

Thickness  of  elecytrolyte 

40  |rm 

20  nm 

Diameter  of  air  feed  tube 

12  mm 

1.4  mm 

Thickness  of  air  feed  tube 

2  mm 

0.1  mm 

Thickness  of  interconnect 

40  |rm 

70  |.im 

perature  distributions  of  a  single  tubular  cell,  in  an  object- 
based  fashion.  The  schematic  image  of  the  framework  is 
shown  in  Fig.  2.  In  this  framework,  models  that  express 
physical  entities  are  constructed  separately  from  models  that 
simulate  physical,  chemical  and  electrochemical  phenom¬ 
ena  occurring  in  the  system.  Physical  entity  models  manage 
the  information  on  dimensions,  microstructure  of  cells 
and  their  components,  and  physical  properties  of  chemical 
substances.  Phenomena  simulation  models  implemented  in 
this  study  are  equivalent  circuit  model,  thermal  model  and 
fluid  model.  Equivalent  circuit  model  simulates  the  current 
flow  and  the  voltage  distribution  in  the  cell  using  the 
structural  information  and  the  physical  property  data  from 
the  physical  entity  models  at  a  given  operating  conditions. 
The  simulated  current  flow  distribution  data  are  input  to  the 
fluid  model  and  the  gas  concentration  distribution  data  from 
the  fluid  model  are  input  to  the  equivalent  circuit  model.  The 
heat  generation  due  to  the  current  flow  through  the  cell  is 
calculated  using  current  flow  and  voltage  distribution  data. 
Heat  generation  data  and  gas  concentration  distribution  data 
are  then  input  to  the  thermal  model.  The  temperature  dis¬ 
tribution  data  from  the  thermal  model  are  input  to  the 
equivalent  circuit  model. 


Fig.  2.  Schematic  image  of  the  model  framework. 
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X-axis 


Fig.  3.  Equivalent  circuit  for  sliced  cell. 


2.3.  Phenomena  simulation  models 

2.3.1.  Equivalent  circuit  model 

To  simulate  the  electrochemical  characteristics  of  the  cell, 
we  used  the  equivalent  circuit  as  shown  in  Figs.  3  and  4.  To 
simulate  the  tubular  cell  characteristics,  the  tubular  cell  was 
divided  into  slices  by  the  planes  perpendicular  to  X-axis,  and 
denoted  by  “sliced  cell”.  The  sliced  cell  is  divided  con¬ 
centrically  into  meshes,  and  they  are  denoted  by  “cell 
mesh”.  A  sliced  cell  was  modeled  by  a  circuit  composed 
of  electromotive  forces  and  cell  impedances,  Z.  The  cell 
characteristics  are  assumed  to  be  symmetric  with  respect  to 
the  broken  line  shown  in  Fig.  3.  The  uniform  electric 
potential  along  the  X-axis  was  assumed  for  the  electrical 
contact  because  the  resistance  of  Ni  felt  is  sufficiently  low. 
Formulation  details  are  described  in  Appendix  A. 

2.3.2.  Fluid  model 

Fluid  model  calculates  flow  rate  and  gas  composition 
distributions  at  the  given  temperature  and  current  flow 


1  I-  :  Electromotive  Force 

Fig.  4.  Equivalent  circuit  for  a  cell  mesh. 


distributions.  For  the  calculation,  fluid  compartment  is 
divided  into  slices  in  the  exactly  the  same  way  with  that 
for  the  cell.  The  slices  in  the  anode  compartment,  cathode 
compartment,  and  air  inlet  tube  compartment  are  denoted  by 
fuel  slice,  oxidant  slice,  and  inlet  oxidant  slice,  respectively. 
The  gas  slice  includes  all  the  slices  mentioned  above.  For 
simplification,  the  gas  concentration  is  assumed  to  be  uni¬ 
form  in  each  gas  slice.  Chemical  species  considered  in  this 
manuscript  are  FI2,  O2  and  H20  for  fuel;  02  and  N2  for 
oxidant.  We  assumed  that  all  the  gases  are  ideal  and  incom¬ 
pressible  fluids.  Faraday’s  law  governs  the  change  in  con¬ 
stituents  due  to  electrochemical  reactions. 

2.3.3.  Thermal  model 

Heat  transfer  occurs  through  thermal  conduction,  con¬ 
vection  and  radiation.  Fig.  5  shows  the  schematic  image  of 
heat  transfer  assumed  in  our  model.  The  air  feed  tube  is 


fljoule  flgen 


Fig.  5.  Schematic  image  of  thermal  transport  phenomena. 
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Fig.  6.  Schematic  algorithm  of  transient  simulation. 


divided  into  slices  by  the  same  plane  dividing  the  tubular  cell 
into  sliced  cell.  In  this  study,  temperature  of  each  sliced  cell 
is  assumed  to  be  uniform.  In  addition,  heat  conduction 
between  sliced  cells  was  neglected  because  it  is  considerably 


small  compared  with  convection  and  heat  conduction 
between  gas  and  solid.  Heat  radiation  was  considered 
between  cell  and  air  feed  tube  but  was  not  considered 
between  cells.  The  heat  exchange  between  solids  and  fluids 
due  to  heat  radiation  was  also  neglected  [15].  Formulation 
details  are  described  in  Appendix  B. 

2.3.4.  Simulation  controller 

The  simulation  controller  controls  both  the  equivalent 
circuit  model  and  the  thermal  model  to  simulate  current  and 
temperature  distributions.  Fig.  6  shows  a  flow  chart  for  the 
simulation  process  of  the  transient  characteristics.  In  the 
transient  phenomena,  the  rates  of  electrical,  chemical  and 
electrochemical  processes  were  assumed  to  be  sufficiently 
fast  compared  with  those  of  thermal  processes.  First  the 
distributions  for  the  steady  state  condition  at  given  operating 
condition  is  derived  after  iterating  the  input  and  output  data 
exchange  between  phenomena  simulation  models.  Then  the 
cell  voltage  is  suddenly  changed  to  simulate  the  transient 
characteristics  of  the  cell.  Next  the  current  and  gas  con¬ 
centration  distributions  at  t  =  At  are  calculated  using  the 
new  cell  voltage  and  temperature  distribution  at  t  =  0.  Using 
the  obtained  current  and  gas  concentration  distributions, 
temperature  distribution  is  subsequently  calculated.  Thus, 
the  transient  characteristics  of  electrical  response  of  the  cell 
are  determined  stepwisely. 

3.  Results  and  discussion 

3.1.  Steady  state  characteristics 

The  current  and  temperature  distributions  at  the  steady 
state  were  simulated  for  both  the  standard  cell  and  the  micro¬ 
tube  cell.  The  operation  conditions  are  listed  in  Table  4.  In 


Fig.  7.  I-V  characteristics. 
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Table  4 

Operation  conditions  for  steady  state  simulation 


Temperature  (K) 

1123 

Pressure  (MPa) 

0.1 

Fuel 

97%  H2,  3%  H20 

Oxidant 

Air  (21%  02,  79%  N2) 

Fuel  Utilization 

0.80 

Air  Utilization 

0.25 

this  study,  temperatures  of  the  inlet  fuel  and  the  oxidant  at 
the  closed  end  of  the  cell  were  assumed  to  be  the  same. 

I-V  characteristics  of  the  standard  cell  are  shown  in  Fig.  7. 
The  temperature  at  the  middle  of  the  cell  is  defined  as  the 
operation  temperature.  I-V  characteristics  reported  by  Sin- 


ghal  [14]  are  also  shown  in  Fig.  7  for  reference.  The 
simulated  I-V  characteristics  for  the  standard  cell  at  890- 
960  °C  show  the  similar  tendency  as  that  of  the  reported 
characteristics  at  900  °C.  However,  the  deviation  between 
experimental  data  and  simulation  data  indicates  that  the 
future  revision  of  the  electrical  and  electrochemical  cell 
characteristics  is  necessary. 

The  current  and  temperature  profiles  along  the  X-axis  for 
the  cells  are  shown  in  Figs.  8  and  9.  The  left  hand  vertical 
axis  is  the  amount  of  current  that  flow  through  the  sliced  cell. 
For  both  cells,  the  temperature  increases  and  then  decreases 
along  the  X-axis.  In  the  micro-tube  cell,  the  current  increase 
along  the  flow  direction  as  temperature  increases,  and 
decreases  before  the  temperature  decreases.  On  the  other 
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Fig.  8.  Current  and  temperature  distributions  at  the  steady  state  for  the  standard  cell. 


Length  [cm] 

Fig.  9.  Current  and  temperature  distributions  at  the  steady  state  for  the  micro-tube  cell. 
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hand,  in  the  standard  cell,  the  current  decreases  along  the  X- 
axis  in  spite  of  temperature  increase.  This  is  because  that  the 
effect  on  Vout  due  to  decrease  of  fuel  concentration  is  larger 
than  that  due  to  the  decrease  of  overpotential  and  ohmic  loss 
as  the  temperature  increases.  The  temperature  difference 
between  the  maximum  temperature  and  minimum  tempera¬ 
ture  is  almost  same  for  both  cells.  Therefore,  the  temperature 
gradient  per  unit  length  is  steeper  in  the  micro-tube  cell 
compared  with  that  in  the  standard  cell.  However,  Kendall 
and  Palin  [3]  demonstrated  that  the  micro-tube  cell  could 
generate  power  without  any  crack  throughout  the  operation 
on  the  severe  condition  that  the  cell  was  directly  heated  by  a 
butane  burner.  These  results  indicate  that  the  stability  of 
micro-tube  cell  is  not  due  to  the  small  temperature  differ¬ 
ence,  but  its  resistance  to  stress.  When  the  temperature 
gradient  in  the  micro-tube  cell  causes  some  problems  never¬ 
theless  of  the  reported  results  [3],  the  extension  of  the  cell 
length  will  mitigate  the  problems. 

The  rates  of  each  type  of  heat  transfer  for  the  sliced  cell  at 
the  center  of  the  standard  cell  are  shown  in  Fig.  10.  Using  the 
obtained  temperature  profile  for  the  standard  cell,  the  rate  of 
heat  conduction  between  sliced  cells  was  calculated.  The 
calculated  rate  of  heat  conduction  is  less  than  1%  of  the  other 
types  of  heat  flow  rates  shown  in  Fig.  10,  which  indicates 
that  neglecting  the  heat  conduction  between  sliced  cells  is 
reasonable  in  the  simulation.  Fig.  10  clearly  shows  that  the 
rate  of  heat  radiation  from  the  cell  to  the  air  feed  tube  is 
considerably  high.  Consequently,  it  is  important  to  take  the 
heat  radiation  between  cells  into  account  in  the  future  work. 

3.2.  Transient  simulation 

The  transient  electrical  response  and  the  temperature 
profile  of  each  component  were  simulated  as  the  terminal 
voltage  drops  from  0.7  to  0.5  V  at  t  =  0.  The  flow  rates  of 
fuel  and  oxidant  were  kept  constant  in  this  simulation.  In  the 
practical  SOFC  system,  fuel  flow  rates  will  be  controlled  to 


keep  the  utilization  ratio  constant  corresponding  to  the  load 
change.  Therefore,  the  inlet  gas  temperatures  will  change 
because  the  inlet  gases  are  preheated  by  using  the  heat  from 
exhaust  gas.  The  temperature  of  the  exhaust  gas  will  vary  as 
the  load  changed,  therefore,  the  temperature  of  fuel  and 
oxidant  after  passing  through  the  heat  exchangers,  will  be 
different  from  that  for  the  rated  operation  condition.  How¬ 
ever,  inlet  temperatures  of  fuel  and  oxidant  were  kept 
constant  for  simplification  in  the  calculation. 

Parameters  for  the  simulation  and  the  conditions  for  both 
cells  before  and  after  the  cell  voltage  changed  at  t  =  0  are 
listed  in  Table  5.  Figs.  1 1  and  12  show  the  transient  power 
profiles  for  the  standard  cell  and  the  micro-tube  cell,  respec¬ 
tively.  The  time  to  reach  the  new  steady  state  for  the  standard 
cell  is  180  s,  while  15  s  for  the  micro-tube  cell. 

The  transient  response  of  the  micro-tube  cell  is  much 
quicker  than  that  of  the  standard  cell  because  the  heat 
capacity  of  the  micro-tube  cell  is  smaller  than  that  of  the 


Table  5 

Initial  and  final  conditions  and  parameters  for  transient  simulation 


Standard  cell 

Micro-tube  cell 

Input  gas 

Fuel  flow  rate  (mol  s_  ) 

2.4  x  1(T3 

3.2  x 

io-4 

Oxidant  flow  rate  (mol  s-1) 

1.5  x  1(T2 

2.0  x 

10~3 

Fuel  inlet  temperature  (K) 

1123 

1123 

Inlet  Oxidant  temperature  (K) 

1040 

1098 

Initial  condition 

Fuel  utilization 

0.46 

0.21 

Air  utilization 

0.17 

0.08 

Current  density  (mA  cm-2) 

208 

722 

Power  (W) 

303 

19 

Final  condition 

Fuel  utilization 

0.80 

0.80 

Air  utilization 

0.30 

0.30 

Current  density  (mA  cm-2) 

358 

2.750 

Power  (W) 

371 

52 
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standard  cell.  The  quick  time  response  is  preferable  for  the 
application  of  SOFC  systems  to  the  power  supply  for 
vehicles  because  it  reduces  the  capacity  of  backup  battery. 


simulation  show  the  quick  response  to  a  load  change  in  the 
micro-tube  cell.  Because  the  quick  response  to  a  load  change 
is  essential  in  transportation  application,  the  micro-tube  cell 
demonstrated  the  possibility  for  transportation  application. 


4.  Conclusions 


A  model  framework  was  developed  in  an  object-based 
fashion  to  simulate  the  transient  characteristics  of  tubular 
SOFC.  Transient  characteristics  of  tubular  cells  with  the  cell 
diameter  of  22  and  2.4  mm  were  simulated.  The  results  of 
steady  state  simulation  show  that  the  temperature  gradient  in 
the  micro-tube  cell  is  six  times  larger  than  that  of  the 
standard  cell.  On  the  other  hand,  the  results  of  the  transient 


Appendix  A 


The  electromotive  force  of  a  cell,  emf,  is  calculated  by  the 
Nernst  equation: 


emf  =  —  In 
4  F 


PQ2c 

P02,a 


(A.1) 
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po2,  R,  T,  and  F  have  usual  meanings.  The  partial  pressure 
of  oxygen  in  anode  was  calculated  assuming  equilibrium 
condition  in  the  gas  phase. 

The  operating  cell  voltage,  Vout,  is  expressed  as  follows: 

Tout  =  emf  -  r/a  -  ric  -  i?IR  (A.2) 

;;a  and  i]c  are  the  anodic  and  cathodic  overpotentials, 
respectively,  and  r/IR  is  the  ohmic  loss  at  the  electrolyte, 
the  electrodes,  and  the  interconnect.  The  operating  cell 
voltage  was  calculated  by  using  the  equivalent  circuit 
model  considering  the  current  distribution  in  the  tubular 
cell. 

Electrode  overpotential  is  consist  of  activation  overpo¬ 
tential,  >/act,  and  concentration  overpotential,  i/conc. 


'/  =  '/act  +  tlconc  (A-3) 

Eq.  (A.4),  given  by  Costamagna  and  Honegger  [16], 
expresses  the  relationship  between  anodic  activation 
overpotential,  ?7act,a,  for  Ni-YSZ  anode  and  current  den¬ 
sity. 


I  =  IQ a 


CXPl  ^RflzcL, 


(A.4) 


r0a  is  the  anodic  exchange  current  density.  Mogensen  has 
proposed  Eq.  (A. 5)  for  the  anodic  exchange  current  density 
[17,18], 


*'o,a  =  ya 


P  h2o 

o,fuel 


exp 


(A. 5) 


The  dependence  of  cathodic  activation  overpotential  for 
LSM  on  current  density  was  formulated  based  on  experi¬ 
mental  results  by  Fukunaga  [19]: 


'7act,c  =  11  x  10  x  exp 


\(l 

1  > 

Ec 

1123  j 

R 

X  Ir. 


(A. 6) 


Ec  is  the  activation  energy  for  cathodic  activation  over¬ 
potential  and  was  assumed  to  be  200  kJ  mol 

The  concentration  overpotential  is  a  potential  drop  caused 
by  the  concentration  difference  in  the  gas  phase.  In  this 
study,  we  assumed  that  concentration  gradient  in  the  gas 
phase  exists  only  in  the  porous  electrodes. 

Diffusion  through  the  porous  material  is  typically 
described  by  either  or  both  of  the  binary  molecular  diffusion 
and  Knudsen  diffusion.  In  this  study,  both  types  of  diffusion 
were  taken  into  account  to  estimate  the  concentration  over¬ 
potentials  [20]. 

The  binary  molecular  diffusion  coefficient  between  gas 
species  i  and  j  in  free  space  is  estimated  by  the  Chapman- 
Enskog  theory  [21]: 


Dij  =  0.0018583 


1  \  */2  7’3/2 

Mj)  P°7jQij 


(A.7) 


We  note  that  the  units  for  total  pressure,  p ,  is  atm  and  a 
is  used  to  denote  the  collision  diameter  only  for  this 


equation.  The  unit  for  collision  diameter  is  A.  Q ,,  is 
the  collision  integral  based  on  the  Lennard-Jones 
potential. 

Eq.  (A. 8)  expresses  the  Knudsen  diffusion  coefficient  for 
gas  species  i. 

Dik  =  97.0  (A.  8) 

V  M-i 

An  effective  diffusion  coefficient  considering  both  binary 
molecular  diffusion  and  Knudsen  diffusion  was  calculated 
as  follows  [20]: 

— —  =  (- —  +  —)  ( (A.9) 
A(e ft)  \Dij  DikJ  \sj 

Tortuous  path  through  the  electrode  microstructure  was 
taken  into  account  by  tortuosity,  £,  and  porosity,  e. 


Appendix  B.  Enthalpy  balance  equations 

The  enthalpy  balance  equation  can  be  described  as  fol¬ 
lows. 

dT 

pCpAV  —  =  ^total  (B.l) 

dr 

p  and  Cp  are  density  and  molar  heat  capacity,  respectively. 
AT  is  the  control  volume  of  each  slice. 

For  the  steady  state,  left  hand  term  in  Eq.  (B.l)  is 
eliminated. 

In  each  gas  slice,  gtotal  can  be  expressed  as  Eq.  (B.2). 

*7 total  —  A  (j  conv  H-  ^<7cond  — g,cond  T  ^^react  (B.2) 

The  net  enthalpy  gain  by  convection,  At/conv,  can  be 
expressed  as  Eq.  (B.3)  using  the  heat  flux  by  convection, 


qcon\i  expressed  as  Eq.  (B.4). 

^^conv  —  ^con  v(^)  ^conv(^  1) 

(B.3) 

^conv(^)  —  ^ 

(B.4) 

m 


Jm  is  the  molar  flow  rate  of  species  m,  and  H  is  the  enthalpy 
expressed  by  Eq.  (B.5). 

fTs(z) 

Hm(T&(z))  =  H°JT0)  +  I  Cp,mdT  (B.5) 

Jt0 

The  net  heat  gain  by  conduction,  A<r/c0ndct]0n,  can  be 
expressed  as 

AtZcond  =  [rg(z  +  1)  +  Tg(z  -  1)  -  2rg(z)]  ^  (B.6) 

S  and  k,  respectively  are  the  area  and  thermal  conducti¬ 
vity. 

The  net  heat  gain  by  heat  conduction  between  solid  and 
gas,  At/s_g  cond,  can  be  expressed  as 

Ags_ g,cond  =  [Ts(z)  -  Tg(z)]Ss-shj,  (B.7) 
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The  net  enthalpy  gain  due  to  mass  transfer  into  the  electrode 
porous  space,  Aqrcuct,  can  be  expressed  as 

Abreact  =  ~  ^ Ho2(Tox(z))  (for  the  oxidant)  (B.8) 

A<7react  =  - A [#h2 (Tfuei (z) ) -HH2o(Tfuel(z))}  (forthefuel) 

(B.9) 

4  is  the  current  density  through  the  electrolyte. 

In  each  solid  slice: 


^/tolal  —  Aqcond  -\-  A  Cf  \/ 1  oss  T  ^^s— g,cond  T  A  q  r  ad  T  Aq  gen 

(B.10) 

A^v.ioss  is  the  heat  generation  due  to  voltage  losses  at  the 
cell. 

The  heat  generation  term  by  entropy  change  through  the 
electrode  reactions,  Aq„en,  can  be  expressed  as 

A<?gen  =  -  ^  A//H,0 (Ts)  -  emf  ie  (B .  1 1 ) 

The  net  heat  gain  by  radiation,  A qrad,  can  be  expressed 
as 

A<7rad  =f^.(T4ccU  -  T4m)  IB.  12) 

Because  the  length  of  the  fluid  channel  along  the  cell  is 
sufficiently  long  compared  with  the  distance  between  cell 
and  air  feed  tube,  only  the  vertical  heat  radiation  between 
cell  and  air  feed  tube  was  only  taken  into  account.  Heat 
radiation  is  assumed  between  the  air  inlet  tube  and  cell. 

The  heat  transfer  coefficient,  h ,  can  be  expressed  by 
Eq.  (B.13)  assuming  a  fully  developed  laminar  flow  at 
constant  wall  temperature,  when  the  Reynolds  number  of 
fluid  is  less  than  2100  [22]. 


“=3.66 

K 


(B.13) 


d  is  the  equivalent  diameter  of  the  gas  channel.  The  thermal 
conductivities  of  gases  are  taken  from  [23].  For  simplifica¬ 
tion,  the  temperature  dependence  of  thermal  conductivity 
was  not  considered. 
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